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Abstract 

X 

We investigate the time-evolutions of angular distributions for B s decays into final states 
that are admixtures of CP-even and CP-odd configurations. A sizable lifetime differ- 
ence between the B s mass eigenstates allows a probe of CP violation in time-dependent 
untagged angular distributions. Interference effects between different final state configu- 
rations of B s — > D* + D*~ , J/ip<p determine the Wolfenstein parameter n from untagged 
data samples, or - if one uses |V^&|/|Vd)| as an additional input - the notoriously diffi- 
cult to measure CKM angle 7. Another determination of 7 is possible by using isospin 
symmetry of strong interactions to relate untagged data samples of B s — > K* + K*~~ and 
B s — > K*° K*°. We note that the untagged angular distribution for B s — > p°<j) provides 
interesting information about electroweak penguins. 
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1 Introduction 



Within the Standard Model |1J one expects || a large mass difference Am = m#— m^ > 
between the physical mixing eigenstates Bf ("heavy") and ("light") of the neutral 
B s meson system leading to very rapid Amt-oscillations in data samples of tagged B s 
decays. In order to measure these oscillations, an excellent vertex resolution system is 
required which is a formidable experimental task. However, in a recent paper |3] it has 
been shown that it may not be necessary to trace these rapid Amt-oscillations in order 
to obtain insights into the fundamental mechanism of CP violation. The point is that 
the time-evolution of untagged non-leptonic B s decays, where one does not distinguish 
between initially present B s and B s mesons, depends only on combinations of the two 
exponents exp(— Tit) and exp(— Tst) and not on the rapid oscillatory Amt-terms. Since 
the width difference Ar = Th — Tl of the B s -sjstem is predicted to be of the order 20% 
of the average B s width [f|], interesting CP- violating effects may show up in untagged 
rates ||. 

In the present paper we restrict ourselves to quasi two body modes B s — > X\ X2 into 
final states that are admixtures of CP-even and CP-odd configurations. The different case 
where the final states are not admixtures of CP eigenstates but can be classified instead 
by their parity eigenvalues is discussed in ||, where we present an analysis of angular 
correlations for B s decays governed by b — > cus quark-level transitions. If both X\ and X2 
carry spin and continue to decay through CP-conserving interactions, valuable information 
can be obtained from the angular distributions of their decay products. Examples for such 
transitions are B s -»■ £> S *+H D tl) D T{-^ D 7l) and B s -> J/V>(-> l + l~)<p(^ K+K~) 
which allow a determination of the Wolfenstein parameter r] || from the time-dependences 
of their untagged angular distributions as we will demonstrate in a later part of this paper. 
Of course, the formalism developed here applies also to final states where the D*^ mesons 
are substituted by higher resonances, such as B s — > D sl (2536) + D s i(2536)~. For many 
detector configurations, such higher resonances may be preferable over D* , because of 
their significant branching fractions into all charged final states and because of additional 
mass-constraints of their daughter resonances. 

If we use the CKM factor 

R b - TTT7T (1) 



with A = sin#c — 0.22 as an additional input, which is constrained by present experi- 
mental data to lie within the ranee R b = 0.36 ± 0.08 @, 1, i, ?7 fixes the angle 7 in the 



usual "non-squashed" unitarity triangle |10| of the CKM matrix |Tl|] through 

sin 7 = ^-. (2) 
Kb 
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Using the isospin symmetry of strong interactions to relate the b — > s QCD penguin 
contributions to B s -> K*+(-^ nK)K*-(^ nK) and B s -> K*°(-^ nK)K^(^ ttK), 
another determination of 7 is possible by measuring the corresponding untagged angular 
distributions. This approach is another highlight of our paper. The formulae describing 
B s — > K* + K*~ apply also to B s — > p°<j) if we make an appropriate replacement of variables 
providing a fertile ground for obtaining information about the physics of electroweak 
penguins. 

This paper is organized as follows: In Section |2| we calculate the time-dependences of 
the observables of the angular distributions for B s decays into final state configurations 
that are admixtures of different CP eigenstates. The general formulae derived in Section [2] 
simplify considerably if the unmixed B s — > X\ X 2 amplitude is dominated by a single 
CKM amplitude. This important special case is the subject of Section |3] and applies to 
an excellent accuracy to the decays B s — > D* + D*~ and B s — > J/ip (j) which are analyzed 
in Section [|. There we demonstrate that untagged data samples of these modes allow a 
determination of the Wolfenstein parameter rj, which fixes the CKM angle 7 if is known. 
In Section ||] we present another method to determine 7 from untagged B s — > K* + K*~ and 
B s — > K*°K*° decays. The formulae derived there are also useful to obtain information 
about electroweak penguins from untagged B s — > p°cj) events. Finally in Section the 
main results of our paper are summarized. 



2 Calculation of the time-evolutions 

A characteristic feature of the angular distributions for the decays B s — > X\ X2 specified 
above is that they depend in general on real or imaginary parts of the following bilinear 
combinations of decay amplitudes: 

A)(t)A f (t). (3) 

Here we have introduced the notation 

A f (t) = A(B s (t)^(X 1 X 2 ) f ) = ((X 1 X 2 ) f \H cS \B s (t)) 
A f ~(t) = A{B s {t) -> {X 1 X 2 ) f ) = ({Xt X 2 ) f \H eS \B s (t)) 

for the transition amplitudes of initially, i.e. at t = 0, present B s mesons decaying into 
the final state configurations / and / of Xi X 2 that are both CP eigenstates satisfying 



(CV)\{Xi X 2 ) f ) = i&KXiXa),) 
(CV) \(X 1 X 2 ) f ) = rjfJiX.X^f) 



(5) 



with ^cp'^cp e { — !)+!}■ Here / and / are lables that define the relative polarizations 
of the two hadrons X\ and X 2 . The tilde is useful for discussing the case where different 
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configurations of X\ X 2 with the same CP eigenvalue are present. To make this point more 
transparent, consider the mode B s — > J ftp <p which has been analyzed in terms of the linear 
polarization amplitudes |L2| A (t), A\\(t) and A±(t) in Whereas A±(t) describes a 



CP-odd final state configuration, both A (t) and A\\(t) correspond to CP-eigenvalue +1, 
i.e. to A f{t) and Aj(t) in our notation with rj{, p = r/^ p = +1. 

The amplitudes describing decays of initially present B s mesons are given by 



A f (t) = A(B 1 (t)^(X 1 X 2 ) f ) = ((X 1 X 2 ) f \H eS \B 1 (t)) 
A~ f {t) = A(B s (t)^(X 1 X 2 ) f ) = ((X 1 X 2 ) I \H cS \B s (t)}. 

Both in these expressions and in (|]) the operator 



(6) 



H cS = H cS (AB = -1) + H cS (AB = +1) (7) 
denotes an appropriate low energy effective Hamiltonian with 

H eS (AB = +1) = H cS (AB = -1)+ (8) 

and 
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j=u,c V *> j=u,c kfc=l k=3 

,(r) 



h cS (ab = -l) = ^ v PQ j = 4 ] { E QicM + E QkCM \ , (9) 



where t> j = V* r Vjb is a CKM factor that is different for b — > d and & — > s transitions cor- 
responding to r = d and r — s, respectively. The four-quark operators Qk can be divided 
into current- current operators (k G {1, 2}), QCD penguin operators (k G {3, . . . , 6}) and 
electroweak penguin operators (k G {7, . . . , 10}), with index r implicit. Note that these 
operators create s and d quarks for r = s and r = d, respectively. The Wilson coefficients 
Cfc(/x) of these operators, where fj, = 0{m^) is a renormalization scale, can be calculated 
in renormalization group improved perturbation theory. The reader is referred to a nice 
recent review 1113 for the details of such calculations. There numerical results for the 



relevant Wilson coefficients are summarized and the four-quark operators Qk are given 
explicitly. 

Applying the well-known formalism describing B s — B s mixing 0, [151 , a straightfor- 
ward calculation yields the following expression for the time- dependence of the bilinear 
combination of decay amplitudes given in (|3]): 

A*f(t) A f (t) = ((X 1 X 2 ) f \H^\B 8 )*((X 1 X 2 ) f \H dS \B 8 ) (10) 



x 



\g + (t)\ 2 + vL C f g+(t) g-(t) + vL g*+(t) + vL vL \g-(t)\- 
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where 



\g±(t)\ 2 

9+{t)g*_{t) 



1 

4 

1 r 

4 



e -r L t + e -r H t ± 2e -rt cos (A m £) 
r L t _ e -r H t _ 2ie ~rt sin ( Amt ) 



(11) 
(12) 



with T = (r^ + T H )/2. The observables and £p which contain essentially all the 
information needed to evaluate the time dependence of fllCf ), are related to hadronic 
matrix elements of the combinations Q? of four-quark operators and Wilson coefficients 
appearing in the low energy effective Hamiltonian (PD through 



E vfl,(X x X 2 ) f \Q?\B s ) 



it 



*(•) 



J=U,C 



(13) 



where 0^ = 2 arg(V^* Vjf,) is the B s — B s mixing phase. In order to evaluate we 
have simply to replace / in ( |13"D by /. Note that we have neglected the extremely small 
CP- violating effects in the B s — B s oscillations in order to derive (|T0|)- (p~3|) [^]. We shall 
come back to ([H|) in a moment. Let us consider the CP-conjugate processes first. The 
expression corresponding to fllCf ) for initially present B s mesons is very similar to that 
equation and can be written as 



A){t) A f (t) = ((X l X 2 ) f ~\H eS \B s )*((X l X 2 ) f \H eS \B s ) 

\g-(t)\ 2 + vL C f 9 m +(t) g-{t) + 7^ 9+(t) 9-(t) + 4 P vL Qtf \9+(t)f 



(14) 



x 



In the general case the tagged angular distribution for a given decay B s (t) — > XiX 2 
can be written as [16 

f(9,v,t,t)=J2bW(t)gW(9,v,i;), (15) 

k 

where we have denoted the angles describing the kinematics of the decay products of X\ 
and X 2 generically by 9, if and ip. Note that we have to deal in general with an arbitrary 
number of such angles. For quasi two body modes B s (t) —>■ XiX 2 into final states that 
are admixtures of CP-even and CP-odd configurations, the observables b^ k \t) describing 
the time-evolution of the angular distribution (pl|) can be expressed in terms of real or 
imaginary parts of bilinear combinations of decay amplitudes having the same structure 
as ([T0|) . The angular distribution for the tagged CP-conjugate decay B s (t) — > X\X 2 on 
the other hand is given by 



f{9^^-t) = Y,b {k \t)g {k) {9^^) 1 



(16) 
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where the observables (t) are related correspondingly to real or imaginary parts of bi- 
linear combinations like (|14]). Since the states Xi X 2 resulting from the B s and B s decays 
are equal, we use the same generic angles 9, if and ip to describe the angular distribu- 
tions of their decay products. Within our formalism the effects of CP transformations 
relating B s (t) — > (XiX 2 )fj and B s (t) — > (XiX 2 )fj are taken into account already by 
the CP eigenvalues 4p and ?7^ p appearing in ([10| ) and (|T^) and do not affect g ( - k \9, if, ip). 
Therefore the same functions g( k \9,(p,ip) are present in ([15]) and (|16|). 

The main focus of this paper are untagged rates, where one does not distinguish be- 
tween initially present B s and B s mesons. Such studies are obviously much more efficient 
from an experimental point of view than tagged analyses. In the distant future it will 
become feasible to collect also tagged B s data samples and to resolve the rapid oscillatory 
Amt-terms. Then Eqs. ( PT0| ) and ( |Hp describing the corresponding observables should 
turn out to be very useful. 

Combining ( [15] ) and ( |T6P we find that the untagged angular distribution takes the form 



[f(9, <p, i/>- 1)} = f(9, <p, ^ t) + f(9, ip, V;t) = E [b (k) (t) + & w (t)] gW(9, <p, (17) 



As we will see in a moment, interesting CP- violating effects show up in this untagged 
rate, if the width difference Ar is sizable. The time-evolution of the relevant observables 
(t) + b^ (t) behaves as the real or imaginary parts of 



A){t)A f {t)\ = A}(t)A f (t)+A)(t)Af(t) = -((X 1 X 2 ) f -\H cS \B s y((X 1 X 2 ) f \H eS \B s ) 



x 



(l+nLvUftf) ( e ~ 



+ e~ 



+ (vL Cj + vL if) (e 



(18) 



In order to calculate this equation, we have combined (|TUj) with (|TJD and have moreover 
taken into account explicitly the time-dependences of (|TT]) and (|l"2"D. We can distinguish 
between the following special cases: 

• /=/: 



\Mt)\ 
x 



= -\((X 1 X 2 ) f \H^B.)\ 2 
(l + K/l 2 ) (e~ rLt + e~ THt ) +2 V f, Re(0) (e" 



(19) 



/ ^ / and r]l F 



A)(t)A f (t)\ = -((X 1 X 2 ) / |iJ cfr | J B s )*((X 1 X 2 ) / |i/ eff |5 s ) 



x 



1 + £/£/) (e- TLt + e" r -') + 4, (Q + 



(20) 
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/ ^ / and ril P = 

^/(*M/(*) 



;((Xl X 2 ) f -|if eff |i? s )*((Xi X 2 )/|f/" e ff|-B s 



(21) 



x 



90) 



-IV.t 



(«?-&)( 



As advertised, the rapidly oscillating Amt-terms cancel in the untagged combinations 
described by ([Tj|). While the time- dependence of fll9D was given in ||, the explicit time- 
dependences of (^0]) and (pi]) have not been given previously They play an important 
role for the untagged angular distribution (|T 



3 Dominance of a single CKM amplitude 



If we look at expression (13), we observe that £y and £j suffer in general from large 
hadronic uncertainties. However, if the unmixed B s — > Xi X 2 amplitude is dominated 
by a single CKM amplitude proportional to a CKM factor , the unknown hadronic 
matrix elements cancel in (|13|) and both and £/ take the simple form 



Q Q < 2 '°- ■ (22) 

where (pf = (axg(yf r Vjb) — a.rg(V t * s Vtb)^J is a CP-violating weak phase consisting of the 
corresponding decay and B s — B s mixing phase. Consequently, in that very important 
special case, ([Tj|) simplifies to 

1 



(r) 



A){t)A f (t)\ = -\((X 1 X 2 ) I \H cS \B s )((X 1 X 2 ) f \H cS \B s )\e l ^ s f^ 



if 



(23) 



x 



-2i<f> 



(r) 



(r) 



(' 



-rv.i 



-Tnt 



where 5f and 5j are CP- conserving strong phases. They are induced through strong final 
state interaction processes and are defined by 



i(f> 



(r) 



-iSj 



(r) 



(24) 
(25) 



((Xi X 2 )f\H cS \B s ) 
((Xi X 2 ) f\H e fi\B s )* 

Note that the structure of (|24| ) and fl2"5|), which is essentially due to the fact that the 
unmixed B s — > Xi X 2 amplitude is dominated by a single weak amplitude, implies that 

_.,(r) -L^M 

the weak phase factors e *^ and e j cancelled each other in (|23| ) and that only the 
strong phases play a role as an overall phase in this equation. We would like to emphasize 
that such a simple behavior is not present in the general case where more than one weak 
amplitude is present. 



The time-evolution of (^3J) depends only on cos ! 
to deal with the following two cases: 



and sin ! 



since we have only 
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vL 



vL- 



Aj(t)A f (t)] = \((X 1 X 2 ) f \H eS \B s )((X 1 X 2 ) f \H eS \B s )\e i ^- s ^ (26) 



x 



cos . 



if 



Vcp Vcp • 



A)(t)A f (t) 

|((X 1 X 2 ) / |i7 cff | J B s )((X 1 X 2 ) / | J f7 eff | J B s )|e^- 5 /^^ p (e"^* - e"^) sin : 



(27) 

<{r) 



Whereas the structure of (|26"D, in particular the cos 2<jyp term, has already been dis- 
cussed for / = / in 0, to the best of our knowledge it has not been pointed out so far 
that untagged data samples of angular distributions for certain non-leptonic B s decays 



allow also a determination of sin 2(pf with the help of (|27j). These sin20y ; terms play an 
important role if the weak phase <jyp is small. The point is that sin 2<jrp is proportional 
to (trf 3 in that case, while cos2(/> < f' ) = 1 + O (^p 2 ^j. Consequently we obtain up to terms 



t(r) 



J 
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V f CP = vL = +1: 



A)(t)A f (t)] = 2\{(X l X 2 ) f \H cS \B s ){{X l X 2 ) f \H^\B s )\e i{& f- & f ) e- TLt (28) 



Vcp Vcp 



-1: 



A}(t)A f (t)] = 2|((X 1 X 2 ) f -|if cfr | J B s )((X 1 X 2 ) / | J ff efr | J B s )|e^- <5 / ) e- r « t (29) 



Vcp 



-Vcp- 



Aj(t)A f (t) 

2t V l P \((X 1 X 2 ) f \H cS \B s }((X 1 X 2 ) f \H cS \B s }\e l{s f- 5 ^ {e~ T ^ - e 



(30) 



r^A jSr) 



We observe that only the mixed combination fl3"0| ) is sensitive, i.e. proportional, to the 

(r) 

small phase (J)j and allows an extraction of this quantity. These considerations have an 
interesting phenomenological application as we will see in the following section. 
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4 The "gold-plated" transitions B s -> D* + D*~ and 
B s ^ J/ip (f> to extract the Wolfenstein parameter 77 



Concerning the dominance of a single CKM amplitude, in analogy to Bd — > J/ijj K$ mea- 
suring sin 2/3 to excellent accuracy [0] (/3 is another angle of the unitarity triangle |TU|). 
the "gold-plated" modes are S s decays caused by 6 — > ccs quark-level transitions. The 
corresponding exclusive modes relevant for our discussion are B s — > D* s + (-^ D~£~j) D*~(^ 
Djj) and B s — > J/i[)(—* </>(— > K + K~). They are dominated to an excellent accu- 
racy by the CKM amplitudes proportional to v)?' = V* s V c b- Therefore the corresponding 
weak phase (f>j defined after ( p2]) is related to elements of the CKM matrix [|nj through 

= [arg(Ofc) - arg(^)]. (31) 

At leading order in the Wolfenstein expansion || this phase vanishes. In order to obtain a 
non- vanishing result, we have to take into account higher order terms in the Wolfenstein 
parameter A = sin^c = 0.22 (for a treatment of such terms see e.g. || |J) yielding |L8|, [H| 

0« = \ 2 r] = 0(0.015). (32) 

Consequently the small weak phase measures simply the CKM parameter r\ [|, [1^, [19|j . 

Another interesting interpretation of ([H]) is the fact that it is related to one angle 
in a rather squashed (and therefore "unpopular") unitarity triangle p0| . Other useful 
expressions for ([H]) can be found in [21]. If we use the CKM factor Rb defined by (|l]) as 



an additional input, rj fixes the notoriously difficult to measure angle 7 of the unitarity 
triangle PH| . That input allows, however, also a determination of 7 (or of the Wolfenstein 
parameter if) from the mixing-induced CP- violating asymmetry arising in Bd — > J/ip K$ 
measuring sin 2(5. Comparing these two results for 7 (or rj), an interesting test whether 
the phases in Bd — Bd and B s — B s mixing are indeed described by the Standard Model 
can be performed. 

The extraction of the weak phase Eq. (32) from B s — > J/ip<j), D* + D*~, etc. is not 
as clean as that of (3 from Bd — > J/ipKs- The reason is that although the contributions 
to the unmixed amplitudes proportional to V* b V us are similarly suppressed in both cases, 
their importance is enhanced by the smallness of <f>^' versus (3 ||22|| . 

Given that is small, we see that (P8|)-(|3T|) apply to an excellent approximation to 
the exclusive channels B s -> D*+(-^ D+j)D*-(-^ Djj) and B s -> J/^(-> 
K + K~), i.e. to X\ X2 G {-D* + _D*~, J/ifj<f)}. Whereas the angular distribution of the latter 
process has been derived in fl3| , a follow-up note |23| not only examines the angular 



distributions for both processes but also discusses an efficient method for determining the 
relevant observables - the moment analysis - and predicts these observables, thereby 
allowing comparisons with future experimental data. 
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The combination ( |30D enters the untagged angular distribution in the form 

lm{[A* f (t)A f (t)}} = 

-2\((X 1 X 2 ) f ~\H eS \B s )((X 1 X 2 ) f \H eS \B s )\cos(8 f - S f ) (e" r ^ - e~ r ^ ^ 



(33) 



C J 



where / G { || , } and / = _L denote linear polarization states p| O]. In order to 
determine the weak phase cj)£> from (|33|) , we have to know both | {{X\ X 2 ) j\H c s\B s ) |, 
|((XiX 2 )/|i?efi|-Bs)| an d the strong phase differences 5f — 5?. Whereas the former quan- 
tities can be determined straightforwardly from 



\A f (t)\\ 

Mt)\ 2 



2\{{X l X 2 ) f \H cK \B s )\ 2 e-^ t 
2\({X 1 X 2 ) i _\H cS \B s )\ 2 e- THt 



(/e{||,o» 



(34) 
(35) 



the latter ones can be obtained by combining the ratio of (|33|) for / 

by 

lm{[AUt)A ± (t)}} 



\((XiX 2 ] 


|||-ffefr|-Bs)| 


cos(<5^ 




\((XiX 2 ] 


\H c g \B S )\ 


cos(<5^ 


-6 ) 



lm{[A* (t)A ± (t)]} 

with the term of the untagged angular distribution corresponding to K 



and / = given 



(36) 



23 



Re{[A* (t)A {] (t)}} = 2|((X 1 X 2 )o|Fe ff |5 s )((X 1 X 2 )|||iJ cff | J B s )|cos(5|| -5 )e~ rit . (37) 

Consequently the angular distributions for the untagged B s — > D* + (— > Dfj) D*~(— > Djj) 
and B s — > J/ , 0(— > l + l~) 0(— > K + K~) modes allow a determination of the weak phase <f>(?\ 

The rather complicated extraction of the strong phase differences 5f — 5j outlined 
above, which is needed to accomplish this task, can, however, be simplified considerably by 
making an additional assumption. In the case of the color-allowed channel B s — ► D* + D*~ 
the factorization hypothesis p5|, |26|, which can be justified to some extent within the 
l/TVc-expansion p7j , predicts rather reliably that the strong phase shifts are mod ir. 
This prediction for the strong phases can be tested experimentally by investigating the 
angular correlations for the 5't/(3)-related modes B u d — > D* s + D* u d . Since B s —> J/ip <fi is 
on the other hand a color-suppressed transition, the validity of the factorization approach 
is very doubtful in this case [p8fl . However, flavor SU(3) symmetry of strong interactions 
is probably a good working assumption and can be used to determine the hadronization 
dynamics of B s — > J/ifxj), in particular the strong phase differences 5f — 5f, from an 
analysis of the S't/(3)-related B — > J/ipK* modes |23|, p4| . These strategies should be 
very helpful to constrain (f>^> with more limited statistics. 

Whereas one expects Th < I"/, and a small value of within the Standard Model, 
that need not to be the case in many scenarios for "New Physics" beyond the Standard 
Model (see e.g. [p9[] ). The untagged data samples described by ( p6|) and ( p7| ) allow then 
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only the extraction of cos 2(jr£' and sin 20^ up to some discrete ambiguities. In particular 
they do not allow the determination of the sign of Ar which could give us hints to physics 
beyond the Standard Model. This feature is simply due to the fact that we cannot decide 
which decay width is Tl and Th, respectively, since we do not know the sign of Ar. 
Using, however, in addition the time-dependences of tagged data samples, sin 2<p^ can 
be extracted and the discrete ambiguities are resolved. With the help of the observables 
corresponding to (^7]) even the sign of Ar can then be extracted, which was missed 
in a recent note [p9| . In general, the ambiguities encountered in studies of untagged 
data samples are resolved by incorporating the additional information available from 
Amt— oscillations. 



5 A determination of 7 using untagged data samples 

of B s K* + K*~ and B s K*°K*B 

After our discussion of some exclusive b — ► ccs transitions and a brief excursion to "New 
Physics" in the previous section let us now consider the b — > uus decay B s — > K* + {— > 
7rK) K*~{-^ ttK) and investigate what can be learned from untagged measurements of its 
angular distribution. Because of the special CKM-structure of the b — ► s penguins ||30| , 
their contributions to B s — > K* + K*~ can be written in the form 

P' f = -\P' f \e i5f p'e iw , (38) 

where / denotes final state configurations of K* + K*~ with CP eigenvalue ?7^ p (see (||)), 
8 p, are CP-conserving strong phases, the CP-violating weak phase has the numerical value 
of 7r and the minus sign is due to our definition of meson states which is similar to the 



conventions applied in 



The penguin contributions include not only penguins with internal top-quark ex- 
changes, but also those with internal up- and charm-quarks [0. Rescattering pro- 
cesses are included by definition in the penguin amplitude P'f. For example, the pro- 



cess B s — > {D* + D*~} — > K* + K*~ (see e.g. [P2| ) is related to penguin topologies with 
charm-quarks running in the loops as can be seen easily by drawing the corresponding 
Feynman diagrams. Although such rescattering processes may affect \P'f\ and 5 P ,, they 
do not modify the weak phase in (j38|). 

On the other hand the contributions of the current- current operators appearing in 
the low energy effective Hamiltonian (0), which are color-allowed in the case of B s — > 
K* + K*~ , have the structure 

T' f = -\T' s \e l8 T' e { \ (39) 
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where 5^, is again a CP-conserving strong phase. Consequently, combining these con- 
siderations, we obtain the following transition matrix element for B s — ► (X\ X 2 ) / with 

X x X 2 = K*+ K* 

(40) 
(41) 



((X 1 X 2 ) f \H eS \B s ) = \P' f \e<> 



where 



\ T 'f\ c usL-s f r ,)^ 



I p'l 
K/l 



Hence the quantitiy £ / defined through ( |13"D is given by 

1 — r/e - * 7 



1 - r f e +i ~*' 



(42) 



Following the plausible hierarchy of decay amplitudes introduced in |3~I|, we expect that 
penguins play - in analogy to B s — ► K + K~ ]33], 0] - the dominant role in B s — > K* + K*~ . 

To evaluate the time-evolution of the observables of the untagged angular distribution 
corresponding to real or imaginary parts of fllSD , we need 1 ± and £~±£j> which are 
given by 



r / ^ r v cos 7 + r / r f 



N 



{'I 



rf) sin 7 



and 



9 + 



iV 



— i 



N 



f,f 

1 — (V~ + r^J COS7 + r*jTf cos 27 

sin 7, 



rjr + r/ — 2rjrj cos 7 



respectively, where 

N fJ ~ _ , , / ^ , ^ , , f 
These combinations of £| and £/ are multiplied in ( |T8| ) by 



l-r%e~ irr -rf e* 1 + r%r f . 



((X 1 X 2 ) f ~\H eS \B s y{{X 1 X 2 ) f \H eS \B s ) = |P;-| iPfle^p'-^')^- 



(43) 
(44) 

(45) 
(46) 

(47) 
(48) 



Here we have used the expression ( flUP to calculate this product of hadronic matrix el- 
ements, which - in contrast to the case where a single CKM amplitude dominates (see 
the cautious remark after fl25|)) - depends also on the weak phase 7 through Nf j. How- 
ever, these factors cancel in ( |I8|) so that we finally arrive at the following set of equations 
describing B s -> (K*+ K*~) f : 
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-1: 



A}(t)A/(t)] = 2|p}| |p;| e i(5 ; 

|l — fr? + 77 J cos 7 + r^rj cos 2 7} 



e Fi * + r-r/ sin 2 7e 



(49) 



^CP 



^CP 



-1: 



A%t)A f {t) 



x 



2|p;-| |p;|e i(5 p'- 5 p' } 



fi /" * , A * 2 1 -Tut 1 * -2 — Frt 

jl — I r; + r/J COS7 + rjTf cos 7 1 e H + r^Vf sm 7 e ^ 



(50) 



^CP 



-^cp 



+ 1: 



z2|p;-| |p;| e i(5 p'- 5 p') 



-r H * 



77 e TLt + r*?rf (e Tht - e rH *jcos7 



sin 7. 



(51) 



The structure of these equations, which are valid exactly, is much more complicated than 
that of (pq)-(pO|) where a single CKM amplitude dominates to an excellent accuracy. Note 
that a measurement of either the e~ Tllt or e~ VLt terms in (49) and (0), respectively, or 
of non- vanishing observables corresponding to flSip would give unambiguous evidence for 
a non- vanishing value of sin 7. 

A determination of 7 is possible if one measures in addition the time-dependent un- 
tagged angular distribution for B s — > K*° K*° which is a pure penguin-induced b — > sdd 



transition. Its time-evolution can be obtained from (49)-(]51"D by setting 



7 = r f 



and 



depends only on the hadronization dynamics of the penguin operators. 

There are two classes of penguin topologies as we have already noted briefly after @: 
QCD and electroweak penguins originating from strong and electroweak interactions, 
respectively. In contrast to naive expectations, the contributions of electroweak penguin 
operators may play an important role in certain non-leptonic P-meson decays because 
of the presence of the heavy top-quark [|35], |3B| (see also [§7]HfHJ). However, in the case 
of the B s — > K*K* transitions considered in this section, these contributions are color- 
suppressed and play only a minor role compared to those of the dominant QCD penguin 
operators. 

If we neglect these electroweak penguin contributions, which has not been done in the 
formulae given above and should be a good approximation in our case, and use furthermore 
the £77(2) isospin symmetry of strong interactions, the B s — > K*° K*° observables can 
be related to the B s — ► K* + K*~ case. In terms of linear polarization states [jP2| , these 
observables fix |Pq|, |P,j|, |Pj_| and cos(8 P , — S p/ ). Since the overall normalizations of 
the untagged B s — > K* + K*~ observables can be determined this way, the e~ VLt and 
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e~ THt pieces of the observables p (t)| 2 ], [|^||0)| 2 ] and Re{[A^(t) A\\(t)}} (see (||)) allow 
another extraction of the CKM angle 7. The remaining observables can be used to resolve 
possible discrete ambiguities. Needless to say, also the quantities r/ and the QCD penguin 
amplitudes Pf are of particular interest since they provide insights into the hadronization 
dynamics of the QCD penguins. A detailed analysis of the decays B s — ► K* + K*~ and 
B s — > K*° K*° is presented in [pH , where also the angular distributions are given explicitly. 

Another interesting application of ([19|) is associated with the decays B s — ► K + K~ and 
i? s — > K° K°. Using again the SU(2) isospin symmetry of strong interactions to relate 
their QCD penguin contributions (electroweak penguin contributions are once more color- 
suppressed and are hence very small), the time-dependent untagged rates for these modes 
evolve as 

\A(t)\ 2 ] = 2\P'\ 2 [(l-2|r|cosp COS7 + |r| 2 cos 2 7)e~ ri * + |r| 2 sin 2 7 e~ THt ] (52) 
and 

~\A(t)\ 2 ] =2\P'\ 2 e~ rLt , (53) 

respectively, where we have used 

r = \r\e ip . (54) 

Here p is a CP-conserving strong phase and \r\ = \T'\/\P'\. In general, there are a lot 
fewer observables in "pseudoscalar-pseudoscalar" cases than in "vector-vector" cases. In 
particular there is no observable corresponding to Re{[Ao(i) A||(t)]}. We therefore need 
some additional input in order to extract 7 from fl5"^) . That is provided by the SU(3) 
flavor symmetry of strong interactions. If we neglect the color- suppressed current- current 
contributions to B + — > 7r + 7r°, which are expected to be suppressed relative to the color- 
allowed contributions by a factor of (9(0.2), this symmetry yields |3l|] 

\T'\ « \ f -^V2\A(B + ^tt+tt )!, (55) 

where A is the Wolfenstein parameter 0, fx and f n are the K- and 7r-meson decay 
constants, respectively, and A(B + — > 7r + 7r°) denotes the appropriately normalized B + — > 
7T + 7r° decay amplitude. Since \P'\ is known from B s — > K° K°, the quantity |r| can be 
estimated with the help of ( |55| ) and allows the extraction of 7 from the part of (|52| ) 
evolving with the exponent e~ Viit . Using in addition the piece evolving with e~ Vht the 
strong phase p can also be determined up to certain discrete ambiguities. Since one 
expects \r\ = 0(0.2) |JT], pi] , it may be difficult to measure the e~ Tllt contribution 



to (p4 ) which is proportional to |r| . The value of 7 and the observable r estimated 
that way could be used as an input to determine electroweak penguin amplitudes by 
measuring in addition the branching ratios BR(5 + — > tt°K + ), BK(B~ — > 7r°i^ _ ) and 
BR(5+ -»■ tt+AT ) = BR(5~ -> n- K^) as has been proposed in ||. 
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Let us finally note that (^9|)-(|5l]) apply also to the mode B s — > p°0, if we perform the 
replacements 



I p> I 
5 f P , 



->'EW\ 



°EWP' 



(56) 



\C' f \ 



I p'EW 



exp 



°EWP< 



where Cf denotes color-suppressed contributions of the current-current operators and 
|Pj EM/ |, &e WP i are related to color- allowed contributions of electroweak penguin opera- 
tors. Similar to the situation arising in B s — > ir°(f), which has been discussed in 



(see also |38|, |39], [40]]), we expect that this decay is dominated by electroweak penguins. 
Consequently its untagged angular distribution may inform us about the physics of the 
corresponding operators. In respect of controlling electroweak penguins in a quantitative 
way by using SU(3) relations among B — > nK decay amplitudes [Q, the CKM angle 7 
is a central input. Therefore the new strategies to extract this angle in a rather clean 
way from untagged B s data samples presented in Sections f| and || are also very helpful to 
accomplish this ambitious task. 



6 Summary 

We have calculated the time-evolutions of angular distributions for B s decays into final 
states that are admixtures of different CP eigenstates. Interestingly, due to the expected 
perceptible B s — B s lifetime difference, the corresponding observables may allow the ex- 
traction of CKM phases even in the untagged case where one does not distinguish between 
initially present B s and B s mesons. As we have demonstrated in this paper, such studies 
of the exclusive b —>■ ccs modes B s — > D* + D*~ and B s — > J/tp <f), which are dominated 
to an excellent approximation by a single CKM amplitude, allow a determination of the 
Wolfenstein parameter 77 thereby fixing the height of the usual unitarity triangle. Using 
the CKM factor Rb oc |T4&|/|K&| as an additional input, 7 can be determined both from r\ 
and from mixing- induced CP- violation in B& — > J/ipKs measuring sin 2(3. A comparison 
of these two results for 7 determined from B s and B d decays, respectively, would allow an 
interesting test whether the corresponding mixing phases are described by the Standard 
Model. 

If we apply the SU(2) isospin symmetry of strong interactions to relate the QCD 
penguin contributions to the b — > uus mode B s — > K* + K*~ and to the b sdd transition 
B s — > K*° K*°, which should play the dominant role there, another extraction of 7 is 
possible from untagged measurements of their angular distributions. Substituting the 
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relevant variables appropriately, the results derived for B s — > K* + K* apply also to 
B s — ► p°0 which is expected to be dominated by electroweak penguin operators. 



We will come back to these decays in separate forthcoming publications |23|, [41| . The 



case of B s decays into final states that are not admixtures of different CP eigenstates 
but only of different parity eigenstates is outlined in ||. There we discuss how angular 
correlations for untagged B s decays governed by b — > cus quark-level transitions allow 
also a determination of the CKM angle 7. 
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